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Abstract 

We calculate the current power spectrum of the gravitational waves created at the big bang (and 
later amplified by the different transitions during the Universe expansion) taking into account 
the present stage of accelerated expansion. Likewise, we determine the power spectrum in a 
hypothetical second dust era that would follow the present one if at some future time the dark 
energy, that supposedly drives the current accelerated expansion, evolved in such a way that it 
became dynamically equivalent to cold dark matter. The calculated power spectrums as well as 
the evolution of the density parameter of the waves may serve to discriminate between different 
phases of expansion and may help ascertain the nature of dark energy. 
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I. INTRODUCTION 



Relic gravitational waves (RGWs, for short) are believed to have their origin in a variety of 
mechanisms at the earliest instants of the big bang. Of particular importance is parametric 
amplification but there are others proposed mechanisms that might also prove very efficients 
in this respect ^,0,0]. Under certain conditions, when the Universe made a transition from 
a stage of expansion dominated by a given energy source to the next, e.g., inflation-radiation 
era, radiation-matter era, these primodial waves experienced amplification. This has been 
widely studied by Grishchuk Q, Allen Q and Maia [6J. Assuming that a mini black hole 
dominated phase existed right after reheating and prior to radiation dominance two other 
transitions may have taken place (namely, inflation-mini black holes era and mini black 
holes-radiation era) with the corresponding amplifications bearing a profound impact on 
the final power spectrum [7[ . Since the cosmic medium is nearly transparent to gravitational 
waves their detection by LISA [l| or some other suitable antenna is expected to yield 
invaluable information on the earliest epochs of the Universe expansion and may allow us 
to reconstruct the whole history of the scale factor. 

Nowadays the observational data of the luminosity of supernovae type la strongly suggests 
that the expansion of the Universe is accelerated at present Q|. In Einstein gravity this 
is commonly associated either to a cosmological constant (vacuum energy) or to a sort of 
energy, the so-called dark energy, that violates the strong energy condition and clusters only 
at the largest accesible scales [Uj . In such a case the present state of the Universe would be 
dominated by dark energy and since it redshifts more slowly with expansion than dust, the 
contribution of the latter is bound to become negligible. Some models, however, propose 
dark energy potentials such that the current acceleration phase would be just transitory and 
sooner or later the expansion would revert to the Einstein-de Sitter law, a(t) oc t 2//3 , thereby 
slowing down (second dust era) |l2|. 



The aim, of this paper is to show how the power spectrum as well as the dimensionless 
density parameter of the gravitational waves created at the big bang, and later enhanced 
by parametric amplifications may help present (and future) observers to ascertain whether 
the expansion phase they are living in is accelarated or not and if accelerated, which law is 
obeying the scale factor. The latter would facilitate enormously to discrimate the nature of 
dark energy between a large variety of proposed models (cosmological constant, quintessence 
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fields, interacting quintessence, tachyon fields, Chaplygin gas, etc) [111] . To this end we 
calculate the power spectrum and energy density of the RGWs when the transitions to the 
dark energy era and second dust era are considered. Obviously, the latter power spectrum 
lies at the future and depending on the model under consideration it may take very long for 
the Universe to enter the second dust era. In general, the scale factor of the FRW metric 
will be of the form a(r]) oc 77™, where 77 is the conformal time and n a constant parameter 
(but different in each expansion phase). Accordingly, we are implicitly assuming that in 
each era there is a constant relation between pressure and energy density. 

Section II briefly describes the different eras of expansion and obtains the Bogoliubov 
coefficients of each transition. Section III calculates the power spectrum of the RGWs at 
the present era of accelerated expansion as well as the power spectrum at the second dust 
era. Section IV determines the evolution of the energy density of gravitational waves within 
the Hubble radius from the present era onwards. Finally, Section V summarizes our findings. 
The notation and conventions of Ref. |^ are assumed throughout. 



II. TRANSITIONS AND COEFFICIENTS OF BOGOLIUBOV 

We consider a spatially flat FRW scenario initially De Sitter, then dominated by radiation, 
followed by a dust dominated era, an accelerated expansion era domiated by dark energy, 
and finally a second dust era. The scale factor in terms of the conformal time reads 



0(77) 
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(—00 < rj < rji < 0), inflation era 



-jj^iv — 2ryi) (rji < rj < 772), radiation era 
*hU iV+ m-2Z l)2 (V2<V<Vs), first dust era (1) 

{ Tit7{m-ti) fa)' (V3<V < Vi), dark energy era 
1 / \ 2 

^ ((I4H4) (77 — 774 + -^jj~ A ) (774 < 77), second dust era 

where I < — 1, rji = 77 + | [(—2// + 1) 7/3 + r] 2 — 4t7i], the subindexes 1, 2, 3, 4 correspond to 
sudden transitions from inflation to radiation era, from radiation to first dust era, from first 
dust era to dark energy era and from the latter to the second dust era, respectively, Hi is 
the Hubble factor at the instant rj = rji. The present time 770 lies in the range [773,774], it is 
to say in the dark energy dominated era. 

Figure ^ schematically shows the evolution of a{rj)H(rj). During the inflationary and dark 
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FIG. 1: Evolution of a{rj)H{rj) in a universe with scale factor given by Eq. We have assumed 
that (14H4 > (12H2. The quantity 77/0 is defined as the instant in the dark energy era in which 
a(r]fo)H(r]fo) = (I2H2, r)f\ as the instant of the second dust era such that 0(77/1)-/? (77/1) = (I2H2, 
and rjf2 as the instant of the second dust era in which a(rjfo)H(r]fo) = 03-^3- 

energy eras a{rj)H(rj) increases with rj, and decreases in the other eras. As a consequence, 
a\H\ results higher than CL3H3. Choosing /, 773 and 774 in such a way that 

x -III / \ 1/2 
04 \ / Gt2 X 



, v - >1, (2) 
0.3/ V a 3/ 

we have that (I4H4 is also higher than a^H^- We assume that 04^4 is lower than a\Hi 
throughout. 

The gravitational waves equation may be interpreted as the massless Klein-Gordon equa- 
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tion 



131 ] . Its solution can be written as 

h v (r Ji x) = J (A {k) hf(k,x)+A\ k) h%\k,xj)d% (3) 



h\j\k.,x) — —=eij(k) - j e tk x , 
J y/iT a(r)) 

»hw+(k 2 -^y ik) (v)=o, (4) 

where A^, A'^ are the annihilation and creation operators, respectively, (k) contains the 
two possible polarizations of the wave and the modes fj,^) (v) satisfy the additional condition 



of the field ^(r?, x) [L3J, [14 |. 



^(k)^*k) ~ ^\k)^[k) = h ( 5 ) 

t 

(k) 



which comes from the commutation relations of the operators A^, A\ k , and the definition 



To study the amplification of gravitational waves at the transition eras we must ascertain 
which of the solutions to (@J) correspond to real particles. An approach to that problem 
is known as the adiabatic vacuum approximation Q|. Basically it assumes that in the 
limit k — > 00 the spacetime becomes asymptotically Mikowskian. In this limit the creation- 
destruction operators of each family of modes will correspond exactly to those associated to 
real particles. 

Following Allen [fj and Maia we shall evaluate the number of RGWs created from 
the initial vacuum state in an expanding universe. Assuming that the scale factor in the 
different eras is given by Eq. the initial state is the vacuum associated with the modes 
of the inflationary stage fii(rj), which are a solution to the equation (j3J) compatible with 
condition (JHJ. Taking into account the shape of the scale factor at this era the modes are 

m = (V^me^k-^x^H^^x), (6) 



r(2) 

— a///, wt an aiuiLiaiv uuii&iaiiu unaoe anu 11 

n 

-3/2 p|. Likewise, the proper modes in the radiation era are 



where jx = krj, ipi an arbitrary constant phase and H_L 2 (x) the Hankel function of order 



fin = ( v ^/2)e*r 1 / 2 x; /2 ffg(x fl ) 1 (7) 

where xr — k(r] — 2t]x) and ip R a constant phase. 

The two families of modes are related by the Bogoliubov transformation 

^(rj) = atiH R {rj) + PiH* R (ri). (8) 



From the continuity of n(rj) at the transition time r/i, it follows that 

ill 
ai ~~ 1 + hf 1 + 2(^)2' A ~ 20^)2' (9) 

where we have neglected an irrelevant phase. The modes with frequency / = 27rk/a(j]i) at 

T]i larger than the characteristic time scale of the transition are exponentially suppressed. 

This time scale is usually identified with the inverse of the Hubble factor at the transition, 

H^ 1 in this case. Therefore, the Bogoliubov coefficients are given by a,\ — 1 and j3\ — for 

RGWs with k > 2ita 1 H 1 , and by Eq. flHJ) when k < 2-na x E x . 

In the first dust era (r] 2 < r] < r] 3 ) the solution to Eq. (J1J) for the modes is found to be 

fjLD = (V^me^k-^x^Hilixn), (10) 

where xjj = k (rj + r\% — 4r]i) and it is related to the radiation modes by 

m(v) = (X2Vd (v) + fon* D (v) ■ ( 1 1 ) 

Consequently, one obtains 



ft= «^ (12) 

for k < 2na{r]2)H2 and a 2 — 1, 02 = for k > 2™ (772)^2, where H 2 is the Hubble factor 
evaluated at r] 2 and r] R2 = V2 — 2r\\. 

Finally, in the dark energy era the solution to equation (|IJ) reads 

M = (V^me^k-^x^H^^), 

where x\ = krji. The Bogoliubov transformation 

Mv) = + PivKv) (is) 

relates the dust and dark energy modes. 
Using the well-known relation 

rW ^ W 2 (-2/)' ,,,,, 



H\%{x) ^e-y-^(-2)V(l + ix), (14) 



valid when a; C 1 and / < —1, it follows that 



«s = [G(-l)- l l a j^)v l Jk l - 2 + 0(k^ 



0: 
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for k < 2710,(773) if 3 and a 3 = 1, (3% = for k > 2ira (773) if 3, where H3 is the Hubble factor 
evaluated at 773 and rj m3 = r/ 3 + 772 — 4r?i. The transition between the first dust era to the 
dark energy era has a time span of the order H^ 1 . This time span is much shorter than 
the period of the waves we are considering and therefore it may be assumed instantaneous 
in the calculation of the coefficients. As a consequence, the time span from this transition 
till today t = to — must be larger than H% , otherwise the transition first dust era-dark 
energy era would be too close to the present time for our formalism to apply. This condition 
places a bound on the value of the redshift — . When I = — 1 the condition is — > 2.72, 
a > 2.25 when I = -2, ^ > 2.15 when / = -3, and > 2 when I -> -00. This bound is 

as '03 ' 123 I 1 

compatible with the accepted values for 22. e [1.5, 11] (see e.g., [la]). Henceforward we will 
consider 22. larger than 2.72 and no larger than 11. 

The solution to equation (J3J) in the second dust era (i.e., the one following the dark energy 
era) is 

where x SD = k(j]-r} A + ■ 

The Bogoliubov coefficients relating the modes of the dark energy era with the modes of 
the second dust era are 

Hi{rj) = a{nsD{v) + ^*sd(v)- 
The continuity of \i at 77 4 implies 

«i = (-i)' l 2l j^r xl if [i 2 + n(i - 2)x u + o(xD] , 

for k < 2-na^H^ and a\ = 1, (3 A = for k > 2-na^H^ where H A is the Hubble factor evaluated 
at the transition time 774. 

III. POWER SPECTRUMS 

In order to relate the modes of the inflationary era to the modes of the first dust era, we 
make use of the total Bogoliubov coefficients ar r and firr- The number of RGWs created 
from the initial vacuum is (N u ) = |/3j> | 2 . Assuming that each RGW has an energy 2hu(r)) 
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(henceforward we shall use ordinary units), the energy density of RGWs with frequencies in 
the range [u(rj) , u(r]) + du)(rj)] can be written as 



dp g {rj) = 2hu(r)) 



2n 2 c 3 



dwirj) 



(N u ) = P{u>{r,))du,(<n), 



(15) 



where P(uj(rj)) = (hw 3 (rj)/ir 2 c 3 ) (N u ) denotes the power spectrum. As the energy density p g 
is a locally defined quantity, it loses its meaning for metric perturbations with wave length 
A = 27rc/uj(ri) larger than the Hubble radius cH~ l (r]). 



A. Current Power spectrum 

To evaluate the present power spectrum one must bear in mind that aoH Q > a^H 3 which 
implies the perturbations created at the transition dust era-dark energy era have a wave 



length larger than the Hubble radius at present 
that a H > a 2 H 2 , it is to say 



17| . One must also consider the possibility 



\ -V / \ 1/2 
a 3 J V a 3/ 

When I = —1 and assuming M ~ 10 4 this condition implies m > 21.5, M > 100 when 

° a 2 S**^ r a 3 ' a 3 

I = —2, and — > 251.2 when / = —3. The values for ^ considered by us are larger than 
2.72 and no lower than 11, consequently we can safely assume aoH < a 2 H 2 . 

For k > 2iraiHi, we find that «T r2 = 1, P? r 2 = 0; in the range 2i^a\H\ > k > 2^02//^, 
the coefficients are ajy2 = «i and Ptt2 = Pi, and finally for k < 2ira(r] 2 )H 2 we obtain 



^2 - ~^T3~2 • ( 16 ) 

8k 3 r]fr] R2 

Thus, the number of RGWs at the present time T)q created from the initial vacuum state 
is (NJ) = \Ptt2? ~ W 6 (^o) for ufo) < 2n(a 2 /a )H 2 , ^~ 4 (r?o) for 2T[(a 1 /a )H l > u(r} ) > 
27i(a 2 /a )H 2 , and zero for u(r} ) > 27i(a 1 /a )Hi, where we have used the present value of 
the frequency, u(r] ) = k/ctQ. 
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In sumary, the current power spectrum of RGWs in this scenario is 
(u(rjo) > 2ic(a 1 /a )H 1 ) , 



n 



(g) 4 Hfu- 1 (2n(a 2 /a )H 2 < u(rj ) < 2n(a 1 /a )H 1 ) , (17) 



^ W (%f (^Y H ^- 3 (2^0 < wfob) < 27r(a 2 /a )H 2 ) . 
While this power spectrum is not at variance with the power spectrum of the conventional 
three-stage scenario (De Sitter inflation, radiation and dust era) 0,0], it evolves differently. 
The power spectrum in the dark energy scenario at rj = 773 formally coincides with Eq. (fTTj) 
but with a 3 substituted by a throughout, and from then up to now waves with 2ira 3 H 3 < k < 
27ia H cease to contribute to the spectrum as soon as their wave length exceeds the Hubble 
radius. By contrast, in the three-stage scenario waves are continuously being added to the 
spectrum. As we shall see in Section IV, this implies that the evolution of the energy density 
of the gravitational waves in the three-stage scenario differs from the scenario in which the 
Universe expansion is dominated by dark energy. 



B. Power spectrum in the second dust era 



e horizon problem posed by an everlasting accelerated 
19^ , models in which the present era is not eternal but 



In an attempt to evade the partic 
expansion to string/M type theories 
it is followed by a decelerated phase as though dominated by cold dark matter (second dust 
era) were proposed 

Here we evaluate the power spectrum at some future time 77 larger than r]f 2 for which 
the waves created at the transition dust era-dark energy era (77 = 773) are considered in the 
spectrum by the first time (see Figure HJ). Let a-7>4 and /?t>4 be the Bogoliubov coefficients 
relating the modes of the inflationary era to the modes of the second dust era. Because of 
condition (J2J) we must consider two possibilities with two different power spectrums. 

(i) If condition (J2J) is not fulfilled, then a\Hi > a 2 H 2 > a^H^ > a 3 H 3 and the power 
spectrum can be obtained from the following total coefficients. In the range k > 2i:a\li\ the 
total coefficients are aiTri = 1 and /3ty4 = 0. For 2na\Hi > k > 2na 2 H 2 , the coefficients are 
«7Y4 = a i an d Prri = Pi, where «i and /?i are defined in Eq. Q. For 2iia 2 H 2 > k > 27ra 4 if 4 , 
o>Tri = «Tr2 and /?Tr4 = Ptv2 where CK7Y2 and /3tv2 are defined in Eq. (fTBj) . 



When > k > 27ra 3 H 3} except for a 3 and p 3} all the cofficients obtained in the 

previous section must be considered when evaluating ar r 4 and Pttl therefore 



a Tr4 = &A a Tr2 ' Pi Ptt2 
PttI = Pi a Tr2 + a A Ptt2 [ 

Finally, for 2ira 3 H 3 > k > 2nH(ri) we get 

-2/)! x 2 



-2i)\ 



32 v ~' ' ~ (-/)! 
32 1 ij (-/)! 



Tr4 
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-/)! 



2/ 



V 2 iVR2Vi4 +2 



ViVR2Vu +2 



,21-5 



J-5 



J-5 



(18) 



VlVR2Vm3Vl4 



(19) 



Accordingly, the power spectrum reads 

(u(rj) > 2^(01/0(77))^) , 



P(„) 



h I 01 \ rr4, ,-1 

4^ l^yJ 



— (£)" (^T)"^" 3 (2vr(a 4 /a(77))^4 < u(rj) < 27r(a 2 /a(7,))ff 2 ) 



(27r(a 2 /o(77))ff 2 < 0,(77) < 2vr(a 1 /a(7 ? ))ii 1 ) 



167T*^C^ V 12/ \ a °/ V a (*?) 



^ V (-01 / \ao/ ^aoj \ao) \a J 

x (4?) 2 ^ 10 Hl°- 2l u 21 - 7 (2ir(a 3 /a( V ))H 3 < 0,(77) < 2ir(a A /a(r,))H 4 ) 



„ 7t ,, ,, . // 9n ,\4 / x 16-81 / \2Z-4 / x2/-4+4/« / x 4Z-8-4/Z 
2 7/t ;4?o4Z-6 / (-2/)! ^ ai] 22. \ 1231 I Ml 

\ (-0! / V a V ^aoj ^aoj ^aoj 

x (4?) 4 ^ 10 <^^" 7 (2vr^(7 ? ) < 0,(77) < 2^(03/0(77))^) • 
(ii) If condition (j2J) is fulfilled, then aiifi > 04^4 > a 2 i/ 2 > 03-^3. As in the previous case, 
in the range k > 2 r na\H\ the total coefficients are «y r4 = 1 and Ptti — 0- For 2Ka\H.\ > k > 
2^a^H^ the coefficients are ar r 4 = oli and Pt t i = Pi- In the range 27ra 4 if 4 > k > 27ra 2 if 2 , 
we obtain 

< 4 = + #a - 4 (-I)- 1 1 2 ^;^; 2 ^ +1 ^- 3 , 

O ( f J- ?7l?7;4 

= + offt ^ <g (-1)" re '"^ 



Again, for 27ra 2 f/ T 2 > k > 2na 3 H 3 , the total coefficients are given by Eq. (J 18)) while for 
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2ira^H 3 > k > 27cH(r]) they obey Eq. f)19j) . The power spectrum in this case is 
{uj{ri) > 2vr(a 1 /a(r/))F 1 ) , 

W (^)) 4 ^i w_1 (2vr(a 4 /a(r ? ))^ 4 < 0,(77) < 27r(ai/afa))2Ti) , 



P(w) ~ < 



9?i 



°(l) 
04 



2/- 



-4/ / \ Z-l / \ 1-3+2/1 / \ 2Z-4-4/Z 
22. I I 2a ) [ 2i 



6-2Z, ,2Z-3 



(27r(a 2 /o(7/))i? 2 < wfa) < 27r(a 4 /a(r ? ))i/ 4 ) , 



(20) 



on ,\2 / \ 16-4Z / \ Z-4 / xZ-4+4/Z / \ 2Z-8-4/Z 
9ft ^2Zo2Z-8 f (~2Z)! \ f ai ^ / 03 ^ (03^ f2i N 



(-0 



"(I 



«o 1 v a °y v a °, 

(2Tr(a 3 /a(i]))H 3 < uj{rj) < 2ir{a 2 /a(rj))H 2 ) 



,, „ /, 9 A|\ 4 / \16-8Z / \22-4/ \2Z-4+4/Z/ n 4Z-8-4/Z 
27?t /4Zq4Z-6 / (-2Z)! \ { g ± \ (a2\ I 2a 1 f 2i 



(t 1 ) 



(-Z)! 
4Z-10 



10—41, ,4Z-7 
1 



(27rfT(77) < 0,(77) < 27r(a 3 /a(r ? )) J f/ 3 ) • 



The power spectrum governed by Eq. (j2()j) is plotted in Fig. |2]for different choices of the 
free parameters I, ^ and as well as the power spectrum assuming the three-stage 
model, i.e., non-accelerated phase and no second dust era. 

The shape of the power spectrum given by ()20j) in the range 2ttH(i]) < u(rf) < 
27T (03/0(77)) if 3 is the same in both cases as in this range all the coefficients are present 
in the evaluation of the total coefficients. It is interesting to see how markedly this spec- 
trum differs from the one arising in the three-stage model (dot-dashed line) at low wave 
lengths. 



1. Topological defects 

Up to now we assumed that at the end of the dark energy era the Universe will evolve 
as if it became dominated by dust once again. Nevertheless if the expansion achieved 
in the accelerated phase were large enough, either cosmic strings, or domain walls, or a 
cosmological constant will take over instead. We will not consider, however, cosmic strings 
(whose equation of state is pk = ~\pk) for, as pointed out by Maia jf|, it seems problematic 
to define an adiabatic vacuum in an era dominated by these topological defects since the 
creation-anihilation operators, A^, A', k >., fail to satisfy the conmutation relations (i. e., 
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FIG. 2: Power spectrum given by Eq. (|2U|) for different values of I. The left panel depicts the 
power spectrum at the instant r/ for which = 10 22 (when Q g (r}) is near unity), with ^ = 11, 
and ^ = 10 6 (to fulfill condition ©). The right panel shows the power spectrum at the instant r\ 
for which ^1 = 10 16 (again when £lg(rj) is near unity), with || = 11, and ^ = 10 6 . For the sake 
of comparison the power spectrum of the three-stage scenario (De Sitter inflation, radiation, dust 
with I = 2) is also shown in both panels (dot-dashed line). Notice the difference in slopes at the 
dust era. 

condition (jJJJ) in the range of frequencies where one should expect RGWs amplification. In 
short, our approach, as it is, does not apply to this case. 

As for domain walls (topological stable defects of second order with equation of state 
Pdw — —\pdw and energy density that varies as a _1 (?7) -see e.g., j^, 21]), once the dark 



energy evolved as pressureless matter at r\ = 774 the scale factor may be approximated by 



a(r) > rj 4 ) = 4a 4 (a 4 if 4 ) ( -77 + r/ 4 + 



2 



-2 



so long as a -1 ^> a~ 3 . That is to say, for 77 > ?] 4 the expansion of the Universe is again 
accelerated whereby a(i])H(r]) resumes growing. The RGWs will be leaving the Hubble 
radius as soon as cH~ l becomes smaller than their wave length, and eventually none of 
them will contribute to the spectrum. 

Finally, we consider the existence of a positive cosmo logical constant A with = A/ (871) 
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and pa = —p\. After the dark energy dynamically mimicked dust the Universe will become 
dominated by a very tiny cosmo logical constant. The corresponding scale factor is 



a(r] > 774) = H A 1 ^-77 + 7/4 + -^Jfj 



Once more, the expansion is accelerated and RGWs will leave the Hubble radius and even- 
tually none of them will contribute to the spectrum. 

IV. ENERGY DENSITY OF THE GRAVITATIONAL WAVES 

By integrating the power spectrum P{oS), defined by Eq. (|15jh one can obtain the energy 
density of the RGWs in terms of the conformal time. Its current value, evaluated from Eq. 
(fTTj) . can be aproximated by [^j 

k / \ 2 / \ 8 

n I a \ 1 di \ TT e, n _ TT N _2 



-32^^; isj fff(2rfW (21) 

To study the evolution of p g {rf) from this point onward we first consider the case > 
a 2 H 2 . In this CcLSG, Pg (77) evolves as (J2TJ) with H{rf) and 0(77) substituted by H and a , 
respectively, till some instant 77/0 in the range 770 < 77/0 < ^4- When 77 ^ 77/0 the RGWs 
with ujijifo) < 2n (02/0(77/0)) H 2 must no longer be considered in evaluating p g as their wave 
length exceed the Hubble radius. Consequently 

h ( d\ \ 4 / Qj\H\ 



< 1 < *> - (^y J J ' (22) 

For 77 = 774, the gravitational waves created at the transition dark energy era-second dust 
era begin contriubuting to p g thereby, 

*<* < , < nn) - *<*) (^y) 4 + ) 2 <*) 

. ; /^y- 1 / 03y- 3+2/ ' /^y^ / y ; - 6 ^ (2^(77))^ 



.<W W W V «4 / 1 -21 + 2 ' 

where Pgirj^) corresponds to Eq.([22|) evaluated at 77 = 774. For 77 = 77/1 > 774 where 77^1 is 
defined by the condition a(r)fi)H(r)fi) = a 2 H 2 (see Figured]), the garvitational waves which 
left the Hubble radius at 77/0 reenter it, therefore 

*< W < ,< W )* ft(w )(^)V ; £^(^) , (|) ,M (2.) 

\ W / x i-4+4/Z / x 2J-8-4/J / / x\2Z-10 / rr/ \\2Z-6 

a 2 \ f«3\ H w-2i( 27TH (W 



ao J \ao J \ a o J V °4 / ' -2Z + 6 
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where p g (?7/i) corresponds to Eq.(j2Sj) evaluated at r\ = rjfx- For 77 = 77/2 > with 77/2 defined by 
0(77/2)^(77/2) — ^3-^3 (see Figure HJ), the gravitational waves created at the transition first 
dust era-dark energy era have wave lengths shorter than the Hubble radius for first time, 
and from that point on the density of gravitational waves can be aproximated by 

»(■> > ^-»^)(^) 1 + S«'- 6 (^) 1 (S) 16 " (25 » 



21-4 / \ 2Z-4+4/Z / \ 41-8-4/1 / , n \ 4Z-10 / t_t/ \\4(-6 

_j f^i) | h 10-4i ( 27tH (V)) 

a J \a J \a J \ a 4 J 1 -4/ + 6 ' 

where p g {j]f 2 ) corresponds to Eq.(j2U) evaluated at r\ = r]f 2 . 

In the case that a^H^ > (I2H2, Eq. (f2~T|) dictates the evolution of p g between 770 till 774. 
Then, from 774 till t]f 2 , p g obeys Eq. (|2"4*j) (note that 7//1 cannot be defined in this case) and 
from rjf 2 onwards p g obeys Eq. (J23j) . 

A natural restriction on p g (r]) is that it must be lower than the total energy density of 
the flat FRW Universe 

where mpi stands for the Planck mass. 

It seems reasonable to consider Hi ~ 10 38 s _1 as it corresponds to the GUT model of 

actor with the scale 
and the current 



inflation |3j, |22j. The redshift ^, relating the present value of the scale 



factor at the transition radiation era-first dust era, may be taken as 10 



value of the Hubble factor H is estimated to be 2.24 x 10 s 231. It then follows 



* = 1.50xlO- 9 ^V" + ", F 2 = 2.24xW^ " V 
a \a 3 J \a 3 

\a 3 J \a J 

— — i8 te) f (s)" 1+f -' 

where we have used the relation 

, 3/2 
<24 x 



H(V > V4) = (-gyj Hi, (26) 

valid in the second dust era (see Eq. ((TJ)) with ajiHjx = a 2 H 2 and aj 2 Hf 2 = a 3 H 3 to 
evaluate if/i (only defined if 04-^4 > a 2 H 2 ) and Hf 2 , respectively. In our model, there are 
only three parameters, namely I < —1, 2.72 < ^ < 11, and 
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We are now in position to evaluate the evolution of the dimensionless density parameter 
£l g (r)) = p g (i])/ p(i]). Its current value is P 



n 



go 



1 Gh f a \ f at . 6 4 
48tt3 c s I aj \a ] ' 



which in our description boils down to 



tt g0 ~ 2.00 x 10~ 13 ( — 



-2+ 



(27) 



Q„o is much lower than unity for any choice of — and I in the above intervals. At later times 



'.yo 

Q g evolves as 



Qg(v > Vo) = fi. 



90 



a 



(28) 



where we have used H(rj4 > rj > r/o) = (^f^J ' H . It is obvious that ^(77) is a 
decreasing function of rj. For rj > rj Q we shall distinguish the two cases mentioned in the 
previous section. 

When condition (j2J) is fulfilled, the evolution of Q g is given by Eq. (|2~£J) till 77 = r]f . 
Then, p g changes in shape from rj = r]f till 77 = 774, as we have seen. Consequently, 



2 Gh fa 



-1+ 



r a 2 / 0(77) 



(29) 



x In 



1 _JL 

4 2i 



ao J 1 ° 



is a decreasing function of 77. Finally, from 774 on, f2 9 evolves as 

a 4 



a{ri) 



24 Gh p l+2l ^2i-t 

7T C 5 



07M\ 2 / \ 8-4Z 



a 



(30) 



x 



, I-l / \ z— 3+2/z , 



a 



«o 



a 



0(77) 



-2/(2 



7T 



-2/ + 2 



(fffa)) 



2Z-4 



Mg{Vf2 > r)>Vfl)- ^g(Vfl) 



0(77) 



+ 



24 Gfr 



2Zq2Z- 



P2 



Z-4 / \ Z— 4+4/i 
(i2 \ / ^3 \ / 04 



a 



7T C J 
2Z— 8— 4/Z 



(-2/)!\V«i xl6 " 4i 



0(77) 
a 1 



(-1)1 

2Z-10 



«0 



(31) 
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(32) 
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0(77) 
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-4/ (27T) 



4/- 



-4/ + 6 
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As follows from ()26|) and (|30 p -()32 p . in each case, the first term redshifts with expansion while 
the second term of Q g grows with rj (recall that I < —1). Conditions Q g (rjfi > i] > 774) <C 1 
and Q g (r]f 2 > r\ > rjfi) <C 1 become just Q, g {r] = 77/1) <C 1 and fl g (r] = r]f 2 ) <C 1, which 
are true in either case whatever the choice of parameters. Finally, from Eq. (J32j) we may 
conclude that at some future time larger than i]f 2 the condition 9 < 1 will no longer be 
fulfilled and the linear aproximation in which our approach rests will cease to apply. 

In the opposite case, when condition (J2J) is not fulfilled, fl g also grows following Eq. (|2"8j) 
till rj = 7/4, in the interval 774 < 77 < 77^2 it grows according to Eq. (|3~Tj). and according to 
Eq. (|H2*|) from 77/2 onwards. Our conclusions of the first case regarding the evolution of Q g 
during these time intervals still hold true. 

The behavior of the density parameter Q„ differs from one scenario to the other. In the 

n 

three stage model (inflation, radiation, dust) fl g remains constant during the dust era |5|]. In 
a four stage model, with a dark energy era right after the conventional dust era, Q g sharply 
decreases in a / dependent way during the fourth (accelerated) era since long wavelenghts 
are continuously exiting the Hubble sphere l^J, see Fig. 01 

By contrast, if a second dust era followed the accelerated (dark energy) era, Q g would grow 
in this second dust era because long wavelenghts would continuously be entering that sphere, 
see Fig. HJ This immediately suggests a criterion to be used by future observers to ascertain 
whether the era they are living in is still our accelerated, dark energy-dominated, era or a 
subsequent non-accelerated era. By measuring Q g at conveniently spaced instants they shall 
be able to tell. Further, if that era were the accelerated one and the Q g measurements were 
accurate enough they will be able to find out the value of the parameter I occurring in the 
expansion law given by Eq. The lower /, the higher the slope of Q g in the second dust 
era. 

One may argue, however, that these observers may know more easily from supernova 
data. Nonetheless, if this epoch lays in the far away future it may well happen that by 
then the ability of galaxies to generate stars (and hence enough supernovae) has seriously 
gone down and as a result this prime method might be unavalaible or severely impared. At 
any rate, even if there were plenty of supernova, the simple gravitational wave method just 
outlined could still play a complementary role. 



16 




-30 -25 -20 -15 -10 -5 -19.2 -1.5 -1 -0.5 
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FIG. 3: Left panel: Evolution of the density parameter £l g with the scale factor from the beginning 
of the radiation era till the present time. The solid line shows the density parameter predicted by 
the three-stage model (De Sitter inflation, radiation, dust). The dashed and dotted lines depict 
the density parameter predicted by a four-stage model having a dark energy era right after the 
dust era with I = —3, —1, respectively, and ao/a^ = 11. Note that the four-stage model predicts 
a lower present value for Q g than the three-stage one as 9 decreases, in a I dependent fashion, 
during the dark energy era. The right panel is a blow up of the region in which the four stage 
models (I = —3, —1), notably differ from the three stage one. 

V. DISCUSSION 

We have studied the power spectrum and the energy density evolution of the relic grav- 
itational waves generated at the big bang by considering the transitions between sucessive 
stages of the Universe expansion. In particular, we have studied the effect of the present 
phase of accelerated expansion as well as a hypothetical second dust phase that may come 
right after the present one. As it turns out, the power spectrum at the current accelerated 
era Eq. (jl7|) formally coincides with the power spectrum of the usual three-stage scenario. 
As a consequence, measurements of P{uj) will not directly tell us if the Universe expansion is 
currently accelerated (as we know from the high redshift supernove data) or non-accelerated. 
However, the density parameter of the gravitational waves evolves differently during these 
two phases: it stays constant in the decelerated one and goes down in the accelerated era 
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FIG. 4: Evolution of the density parameter Q g given by Eqs. (|27|) - (|32|) with the scale factor 
in a five-stage scenario (De Sitter inflation-radiation-dust-dark energy-second dust era) from the 
beginning of the dark energy era. The second dust era starts at the time in which the graphs 
attain their minimum. The solid, dotted and dashed lines correspond to I = —1, I = —1.5 and 
I = —2, respectively. The density parameter decreases in the dark energy era and increases in 
the second dust era till its value is comparable to unity. From that point on our approach (being 
linear) ceases to apply. We have chosen a dark energy era which begins when 03/00 = 1/11 and 
ends when 04/00 = 10 6 . 

in a I dependent manner. Therefore, the present value of fl g may not only confirm us the 
current acceleration but also may help determine the value the parameter / -see Fig. |3} and 
hence give us invaluable information about the nature of dark energy. 

In the far away future measurements of P{uj), if sufficiently accurate, will be able to tell 
if the Universe is still under accelerated expansion (driven by dark energy) or has entered a 
hypothetical decelerated phase (second dust era) suggested by different authors [l^. This 
may also be ascertained by measuring the density parameter of the gravitational waves at 
different instants to see whether it decreases or increases. 
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